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We propose a scenario of brane cosmology in which the Peccei-Quinn field plays the role of 
the infiaton and solves simultaneously many cosmological and phenomenological issues such as the 
generation of a heavy Majorana mass for the right-handed neutrinos needed for seesaw mechanism, 

MSSM /r-parameter, the right amount of baryon number asymmetry and dark matter relic density at 
the present universe, together with an axion solution to the strong CP problem without the domain 
wall obstacle. Interestingly, the scales of the soft SUSY-breaking mass parameter and that of the 
breaking of D(1)pq symmetry are lower bounded at 0(10) TeV and 0(10^^) GeV, respectively. 


INTRODUCTION 

Probably, the most compelling source for the density 
perturbations in the present universe would be the classi- 
calized quantum fluctuations of the infiaton [1, 2]. How¬ 
ever, in addressing the spectral features of the density 
perturbations observed by Planck satellite [3], the flat¬ 
ness problem (called 77 -problem) of the infiaton potential 
[4, 5] turns out to be a big obstacle in the attempts of 
a simple realization of inflation based on an UV theory, 
e.g., supergravity or string theories with Einstein gravity. 

Actually, the 77 -problem can be removed rather easily if 
the infiaton path is a simple compactified trajectory in a 
multi-dimensional field space without resorting to trans- 
Planckian excursions for the infiaton field [6-15] (see also 
[16]). Notably, the resulting inflation is effectively the 
same as single-field slow-roll inflation. 

Another interesting possibility for tackling the 77 - 
problem is to introduce an extra-dimension [17]. Depend¬ 
ing on the relative size of the energy density on the brane 
as compared with the brane tension, the expansion rate 
on the brane can be much larger than the one expected 
in Einstein gravity. This can resolve the 77 -problem, al¬ 
lowing inflation even with a steep potential [18]. 

In this letter, we show that the Peccei-Quinn field [19] 
of the DFSZ axion model [20, 21] in the supersymmetric 
framework is a natural candidate for the infiaton which 
can also trigger some necessary post-inflation cosmology 
such as baryogenesis and dark matter,when considering 
brane cosmology. 


INFLATION ALONG A FLAT DIRECTION ON A 
BRANE 

In Randall-Sundrum type II brane world scenario [17], 
after the bulk contribution to the brane is diluted away, 
Friedman equation on the brane can be written as [22] 

3H^M^ =p{l + £) ( 1 ) 

where Randall-Sundrum condition for vanishing cosmo¬ 
logical constant is assumed, H is the expansion rate on 


the brane, Mp = 2.4 x 10^® GeV is the reduced Planck 
mass, and p is the energy density on the brane other than 
the brane tension A/2 which is constrained to be larger 
than about (2.3 TeV)^ from the validity of the (3 -I- 1)- 
dimensional Newtonian gravity on scales larger than 
about 44/im [23, 24]. It implies that, if V = p/K ^ 1 
with p dominated by the potential energy V of a scalar 
field, a long enough inflation is possible even with a steep 
potential which would not allow or sustain an inflation¬ 
ary epoch in Einstein gravity. The slow-roll parameters 
of such an inflation on the brane are given by [25] 
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where eE,??E and are the conventional slow-roll pa¬ 
rameters in Einstein gravity. The spectral index and its 
running can be expressed as 
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Notice that the spectral index is given by the same ex¬ 
pression as in the case of Einstein gravity. The power 
spectrum and tensor-to-scalar ratio yield [25] 


P'R.= {^ + Pn,E, rr = ^ (6) 

with Pn,E = V/Mp€p), while the number of e- 
foldings is 



with (j) being the infiaton field. For the horizon scale 
of the present universe ~ 6000Mpc, the e-foldings 
expanded at the time the Planck pivot scale (fc* = 
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0.05Mpc exits the horizon during inflation becomes 
[26] 
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where So is the total entropy of the present observable 
universe, the subscript 4 ’ in the right-hand side denotes 
a value at the end of inflation, and Tj is the decay tem¬ 
perature of the inflaton. Note that the last term in the 
second line of Eq. ( 8 ) is the new contribution coming 
from the enhancement of the expansion rate. Also, in 
order not to over-produce KK-gravitions, < (A/2)^/^ 
is demanded [27]. 

Generically, thanks to the non-renormalization theo¬ 
rems, the potential E of a supersymmetric flat direction 
($ = 0 e®®/-\/ 2 ) is dominated by a mass term until the 
field gets close to the true vacuum or before it is lifted 
up by higher order terms. For a symmetry-breaking flat 
direction, the potential can be written as 


The power spectrum can be expressed as 

^ 12^2|^|3 (^) (15) 

and P-R = 2.142 x 10“® which with 77 = —1.7 x 10“^ [3] 

leads to 

m/ 4 o ^ 10 “^ (16) 

Also, from Eqs. (14) and (16), the scale of brane tension 

is found to be 

Since A is constrained as (A/2)^/^ > 2.3TeV, we hnd 

e>(10)TeV<m<ms (18) 

with TTis being the scale of the soft SUSY-breaking mass 
parameter, and 


((>0 > 0(10^1) GeV (19) 

It is very interesting to notice that, the Peccei-Quinn field 
is a natural and compelling candidate of (f). 
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THE MODEL 


where Vq is set for vanishing cosmological constant at 
the true vacuum, m?(> 0 ) is a mass-square parame¬ 
ter, and ‘... ’ represents higher order terms suppressed 
by a very large mass scale, e.g., Planck scale. Roughly, 
Vo with (j)o being the vacuum expectation value 

(VEV). Then, for ^ <C (/)o one finds 
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Hence, if ( 7 !) <C 4o for a relevant cosmological scale during 
inflation, the slow-roll parameters of brane inflation are 
such that e <C 77 resulting in 
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Combining with the number of e-foldings in Eq. (7), we 
find 
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Eq. (12) implies 
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Motivated by the observation in the previous section, 
we consider the Peccei-Quinn field as the inflaton. For a 
specific realization, we consider a model defined by the 
following superpotential where gauge group and family 
indices are omitted: 

W = YuQHuU + YaQHdd + YeLHae 
+^XN^ + X^LHuN 

, XxyX^Y 

Mp Mp ^ 

here U{1) Peccei-Quinn symmetry is assumed with 
charges assigned as qpQ{X,Y) = (1,-3). The first line 
of Eq. (20) is the supersymmetric realization of the stan¬ 
dard model (SM) Yukawa couplings. The second line 
triggers the seesaw mechanism when X develops a non¬ 
zero vacuum expectation value (VEV). That is, right- 
handed neutrinos (RHNs) get masses given by 

rriN = Xn{X) (21) 


cj), ~ _ Q 37 ^^ ^ O(0.1)4o (13) 

where for a numerical estimation we used A4_* = 50 and 
77 Ri —1.7 X 10“^ [3] leading to 

V (14) 


where the family index was suppressed. As mx H, 
RHNs are integrated out, replacing the second line of 
Eq. (20) by 
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which generates Majora masses of active neutrinos at low 
energy given by 
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with Vu being the VEV of The third line of Eq. (20) 
is nothing but a supersymmetric realization of the DFSZ 
axion model [28]. The first term of the third line repro¬ 
duces the ^-term of the minimal supersymmetric stan¬ 
dard model (MSSM) as X develops a non-zero VEV, and 
the second term stabilizes X and Y. 

Eor a simple illustration of inflation and post-inflation 
cosmology in the model of Eq. (20), we set X = Y = ^ 
and consider a superpotential, 


W = 


A 

4M^ 


(24) 


where A = 0(0.1 — 1) is a dimensionless constant. Then, 
including the soft SUSY-breaking terms with a negative 
mass-square parameter assumed [36], the scalar potential 
of ($*) becomes 


V = Vn-m^l<Pl^ - — 
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where Vq is again set to get a vanishing cosmological con¬ 
stant at true vacuum, and to^(> 0) and A are soft SUSY- 
breaking parameters. The vacuum expectation value of 
4) is found to be 
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The masses of radial and angular modes at true vacuum 
are 
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COSMOLOGY 

The model of Eq. (20) can accommodate inflation in a 
brane world scenario, and trigger the realization of a suc¬ 
cessful baryogenesis and supply dark matter production 
as follows. 


Inflation 


Although inflation in the model of Eq. (20) takes place 
along X, all its main features can be captured from the 
simplified example of Eq. (24). From our numerical anal¬ 
ysis, we found that, as an example, the following set of 
parameters gives a perfect fit to the Planck data: 

TO ~ 1.4 X 10® GeV, A ~ 0.61, 2V = 10^® ®^ (30) 

which results in 


TOpQ,r = 1.96 X 10® GeV, (j)o = 9.15 x 10“ GeV (31) 

and (A/2)^ ~ 51 TeV. One can also lower (j>o while keep¬ 
ing m/(j)o nearly invariant for a correct value of Pjz. The 
decay rate of the flaton is assumed to be [32] 




Vi ^ 

47r V rn\J \mpQ^J <(>o 


(32) 


where B is the soft SUSY-breaking parameter of the 
Higgs bilinear term in MSSM and iriA is the mass of 
CP-odd Higgs. For niA = 1.03H and /r = to, one finds 
Td ~ 550 GeV leading to A^e.* — 55. Inflation ends at 
(pe ~ 0.8(/)o, and inflationary observables turn out to be 


10®Hr ^ 2.136, ns ~ 0.962, = -6.1 x 10"^ (33) 

am fc 


The tensor-to-scalar ratio is miserably small, undistin- 
guishable zero from an experimental point of view. 

Even if 4> evolves in a 2-dimensional field space, for 
4’ ^ 4’e, the angular direction of $ is much flatter 
than radial direction. Hence, unless the trajectory of the 
inflaton is around the ridge of an angular potential (e.g., 
the potential involving A-parameter in Eq. (25)), the 
non-Gaussianities are expected to be /nl ^ 1, since per¬ 
turbations are not enough to produce significant changes 
on the inflaton trajectories which nicely converge to the 
true vacuum. This is expected to be true as well for 
the real model of Eq. (20) where X plays the role of the 
slow-roll inflaton. Note that Y develops its VEV via the 
A-term of the second term in Eq. (20). So, Y is expected 
to follow its time-dependent vacuum position during the 
whole period of inflation, and therefore it would not af¬ 
fect the adiabatic density perturbations caused mainly 
by the perturbations of X. 

A remark on the initial condition of the inflation along 
PQ-field is in order. If the PQ-field were the only light 
field, we can simply assume a radiation dominant uni¬ 
verse before the onset of inflation. However, there are 
at least two possible candidates of light fields too: GUT 
Higgs and Planckian moduli. For A^/'^ = 0(10®“®) GeV 
which is the case for the parameter set we have chosen for 
illustration purposes, they can also lead inflation while 
they roll out/in from/to the origin. As a simple possibil¬ 
ity to realize PQ-inflation as the last inflation that may 
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be necessary to have a proper post-inflation cosmology 
(e.g., baryogenesis without unwanted relic problem), we 
may assume that the mass scales of those field are slightly 
higher than that of the PQ-field. In this case, the expan¬ 
sion rate at the end of the pre-existing inflation (i/pre) 
could be a little bit larger than the one from PQ-inflation. 
The PQ-field would barely move from the origin once it 
is there. The potential danger of being drifted away by 
quantum fluctuations can be removed if PQ-field couples 
to a light degree of freedom around the origin. Specif¬ 
ically, in the presence of a coupling W D the 

effective mass-square of the PQ-field around the origin is 

( 34 , 

since RHNs would have fluctuations of order Hp,-e/{2Tr). 
Hence, if 

Hpre > 271171/Xn (35) 

the PQ-field can be safely held at the origin. In other 
words, the mass scale of GUT Higgs or Planckian mod¬ 
uli is required to be larger than that of PQ-field by a 
factor about 2-k/Xn. If GUT Higgs causes pre-inflation, 
the reheating temperature is likely to be low, since the 
interaction to light fields will be suppressed by GUT 
scale. In this case, the PQ-inflaton may start rolling 
out soon after the pre-inflation if Hp^e ^ 277771 /X^. 
On the other hand, Planckian moduli may end up at 
a symmetry-enhanced point leading to a rapid reheat¬ 
ing and resulting in a rather high reheating temperature 
(e.g., ^ (6iJp;.gMpA)^^®). In this case, the PQ-infiaton 
could be held at the origin until the temperature of ther¬ 
mal bath becomes comparable to m (a phase of thermal 
inflation [29, 30]). 

Our choice of parameters results in H ^ 
0(10) (A/2)^^^ during inflation. In this case, the 
higher order kinetic terms of the inflaton, caused by 
brane fluctuations, are not suppressed enough, making 
the quadratic kinetic term a rather poor approximation 
[37]. It is non-trivial to properly handle this issue given 
the tachyonic nature of our inflaton, and out of the 
scope of this paper. However, at least, once the brane 
tension is given, enforcing the validity of considering 
only the quadratic kinetic term of the inflaton as we 
did in this work may impose a constraint on the energy 
density of the inflation Vq such that Hubble fluctuation 
is at most of the same order as the brane tension. As 
a simple solution to this issue in our scenario, one can 
increase the size of A-parameter relative to the mass 
parameter m by a factor ~ 0(10) in Eq. (25). In this 
case, the contribution of the A-term to the potential 
for the cosmological scales of interest becomes sizable, 
requiring a smaller in order to obtain the right 
amount of e-foldings. This means that e, becomes 
smaller. As a result, for a given brane tension one can 


lower the expansion rate during inflation by taking a 
slightly smaller mass scale for both parameters, m and 
A in Eq. (25), while keeping inflationary observables 
nearly unchanged. Baryogenesis and the production of 
dark matter discussed in the next subsections are barely 
affected although there will be minor changes in the 
working parameter space. So, we do not pursue the 
details of this case here. 


Baryogenesis 

The dynamical generation of /i can sustain a late time 
Aflleck-Dine leptogenesis, as discussed in Refs. [31-33], if 

= (™i + /2 < 0: -t > 0 (36) 

at scales below an intermediate scale, where 7n\ and 
are respectively the soft SUSY-breaking mass-square pa¬ 
rameters of L and Hu fields. In addition to Eq. (36), we 
also assume 

(37) 

with —being understood as the soft SUSY-breaking 
mass-square of X so that the PQ-field rolls away from 
the origin following the Li7„-flat direction. Note that, in 
the model of Eq. (20) when X is held around the origin, 
the LHu direction is stabilized by a quartic term. How¬ 
ever, as X rolls out and becomes large, the renormaliz- 
able terms in the second line of Eq. (20) are replaced by 
the one of Eq. (22). The inverse of the time scale of the 
evolution of (£), the time-dependent VEV of the LHu 
flat-direction, is 

{£)/{£) = X/X = -r^xH (38) 

where rjx is the slow-roll parameter of the inflaton {X). 
It is much smaller than the mass scale of LHu around 
the potential minimum. Hence, the LHu flat-direction 
would trace its potential minimum as X evolves. 

The lepton number asymmetry at the onset of the an¬ 
gular motion of the AD field is estimated to be 

riL ~ sin(40, + A0) (39) 

TOpQ.r AatAq 

with TOpQ^i- being the mass of the radial mode of X 
around the true vauum, and £e is the VEV of LHu flat 
direction at the end of inflation. Because the oscilla¬ 
tion of X causes time-dependent variations of the mass 
and of £e of LHu, the initial lepton number asymme¬ 
try is suppressed by a factor which we denote as i5s. 
Passing through the electroweak symmetry breaking, the 
asymmetry in the lepton number is converted to the one 
in baryon number via anomalous electroweak processes. 
Then, taking ns Ui/3, the late time baryon number 
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asymmetry after the inflaton decay can be expressed as 

^ j Td f A, V f lo y 

S \mpQ^rJ \XoscJ 

x/(-4i-) ( sin(46>^ + AO) (40) 


where we use 
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where = AatATo and Ap = A^, Xqsc = ckXq is the 
initial oscillation amplitude of X at the end of inflation 
with a = 0(0.1), /(Ap) = X^Iq/ and A0 is 

the relative CP-violating phase between the two lepton- 
number violating terms [33] and assumed to be of order 
one. The precise estimation of Sg requires a heavy nu¬ 
merical siirrulation which is out of the scope of this let¬ 
ter. However, froirr the simulatioir results of Ref. [32], 
we expect 6 s = 0 ( 10 “^ — 10 “^), taking into account 
the larger expansion rate at the time of AD leptogen- 
esis and the weaker preheating due to the smallness of 
the ^-term-induced mass variation of the LH^ direction. 
For the parameter set used in the previous section with 
^ obtain £o ^ 10 ®“^°GeV and 

the factors appearing in the second line of Eq. (40) are 
0(0.1) (once taken together). Hence, the right amount 
of baryon number asymmetry would always be obtained 
if 


suppression of Sir^a- Observations constrain the baryon 
isocurvature perturbations to be /3iso = Ps/ {Pn + Ps) ^ 
0.038 with Pg being the power spectrum of the isocurva¬ 
ture perturbations [34]. This can be interpreted as 

^ ^ 5.34 X 10-5 (45) 

\1 Piso / 

where Hdm and H b are the densities of cold dark matter 
and baryons, respectively. For A^ = ypppj = m as an 

example, Sh-y = 2.37 x {m/I q/\Q~^). Therefore, we 

are safe as regards of baryon isocurvature perturbations. 


Dark matter 

If the inflaton decays well after the freeze-out of the 
lightest neutralino which is likely to be our case, in order 
to avoid dark matter over-production, the inflaton de¬ 
cays into flatinos and neutralinos should be kinematically 
forbidden (for example, see Ref. [35] for the parameter 
space). Hence, dark matter in our scenario is likely to be 
either neutralinos from thermal bath or axions, depend¬ 
ing on the choice of parameters and the initial conditions 
of inflation. 

The energy density stored in radiation after inflation 
until the energy density of PQ-field (ppq) becomes com¬ 
parable to the brane tension is expected to be 

Pr ' Tlh}'^ (46) 


Td = 0(10-5 - 10-2)m (42) 


Notice that, unless 6 s is too small, the right amount of 
baryon number asymmetry can be obtained by adjusting 
(mostly lowering down) Td (or B/tua - see Eq. (32)). 

Since the AD field is a flat direction during inflation, 
there is a potential danger of generating too large baryon 
isocurvature perturbations {Sb~f = <5 In (ns/s)), but they 
are naturally suppressed in our scenario for reasons which 
will be explained below. 

Setting {9i) = 0 without loss of generality, from 
Eq. (39) one finds 

where 6 ir,a is the radial/angular perturbations of LHu 
and we used 69(^ = 6 £al£o- Denoting the rj of Ti7„ around 
the true vacuum as 77 ^’g for the radial and angular modes 
of LHu, the perturbations of each mode at the end of 
inflation can be written as 


6£r 


6e 




(44) 


where 6^1 ^ = H/ (27r). Note that can be larger than 
\rix\ during inflation by a factor of a few, providing a 


where g*(rd) is the number of relativistic degrees of free¬ 
dom at inflaton decay. Hence, the background tempera¬ 
ture during this epoch (Tx) is nearly constant irrespec¬ 
tive of iJ, and it may or may not be lower than the typical 
freeze-out temperature of neutralinos in Einstein gravity, 
Tfz ^ m ^/20 with being the mass of neutralino. 

If Tx > Tfz, for Td < Tfz the abundance of neutralinos 
would be diluted by the entropy production generated in 
the decay of the inflaton. The dilution factor is given by 



32/3 

V 3* (2d) / V 2d / 


(47) 


where Td is the would-be background temperature when 
H ^ Fd in case of no entropy production. The relic 
density of neutralino can then be expressed as 

^DM ( my Y I „io/ my Y ( 

"" Vl02GeV/ VlO^GeV/ \my^J 

(48) 

where we used g*(Tfz) = g*(Td) and Tfz = rriy.(20 for A. 
Hence, for rriy^ ~ nis ~ O{10‘^~^) GeV the abundance of 
neutralinos can match the observed relic density of dark 
matter if T^/my_ ^ 0(10-5 — IQ-^) which may also be 
compatible with that required for baryogenesis. 
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If Tx < Tfz, neutralinos would never be in thermal 
equilibrium and their abundance will be negligible, its 
exact size depending on Tx /Tf^. In this case, dark matter 
would have to be mainly axions. 

As PQ-symmetry is broken during inflation, the axion 
domain wall problem is absent, and the contribution to 
dark matter comes only from misalignment. In this case, 
the axion coupling constant can be much larger than the 
typically allowed upper-bound, since the misalignment 
angle is given as the initial condition of inflation. 

CONCLUSIONS 

In this letter, we proposed an inflation scenario along a 
supersymmetric flat direction in the context of Randall- 
Sundrum type-II brane world. Showing that, the mass 
scale (rn) of the symmetry-breaking supersymmetric flat- 
direction is directly connected to the symmetry-breaking 
scale {(j)o) as m/cj)^ ^ 10“”^, we identified the Peccei- 
Quinn field associated with the axion solution to the 
strong CP-prblem as the natural candidate to play the 
role of the inflaton. We also showed that a successful 
post-inflation cosmology including small active neutrino 
masses, baryogenesis and dark matter can be triggered 
by the non-zero vacuum expectation value of the PQ- 
inflaton via the generation of the heavy Majorna mass 
of the seesaw mechanism and the dynamical generation 
of the MSSM /r-parameter at the end of inflation. Either 
neutralinos or axions can be the main component of dark 
matter without an axion domain wall problem. 

Our scenario does not have ad-hoc parameters ex¬ 
cept the brane tension (A). The fundamental Planck 
scale in five-dimensional spacetime is expected to be of 
0(10®) GeV with a brane tension of 0(10'*“®) GeV. The 
brane tension imposes lower bounds on the scales of soft 
SUSY-breaking mass and breaking of C/(1)pq symmetry 
at 0(10) TeV and 0(10**) GeV, respectively. This sce¬ 
nario can be excluded if primordial tensor perturbations 
are observed, a TeV scale SUSY (meaning A*/* < 1 TeV) 
is found, or the axion coupling constant turns out to be 
lower than 0(10*®) GeV. 

In summary our scenario is a full description of the 
early history of the Universe from inflation to Big-Bang 
nucleosynthesis that may be easily falsifiable at future 
collider experiments. 
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